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ABSTRACT: In the present paper we derive six and seven loop formulas for the anomalous 
dimension of the Konishi operator in N=4 SYM from string theory using the technique of 
Luscher corrections. We derive analytically the integrand using the worldsheet S-matrix 
and evaluate the resulting integral and infinite sum using a combination of high precision 
numerical integration and asymptotic expansion. We use this high precision numerical result 
to fit the integer coefficients of zeta values in the final analytical answer. The presented six 
and seven loop results can be used as a cross-check with FiNLIE on the string theory side, 
or with direct gauge theory computations. The seven loop level is the theoretical limit of 
this Luscher approach as at eight loops double-wrapping corrections will appear. 
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1 Introduction 



The Konishi operator is the simplest nonprotected local operator in J\f = 4 Super- Yang- 
Mills theory. For this reason, its properties serve as a theoretical laboratory in our quest 
for the complete solution of J\f = 4 SYM at any coupling. The relative simplicity of its 
structure and properties allows us to probe them using various methods - including direct 
gauge theoretical computations and computations on the string side of the AdS/CFT cor- 
respondence. On both sides of the AdS/CFT duality there are various ways (sometimes 
conceptually quite different) of performing these computations and the agreement of the 
computations was always a strong confirmation of the respective methods. Apart from 
the above mentioned comparisons between different computations for the Konishi opera- 
tor, there are also sometimes some quite nontrivial internal consistency checks within one 
particular method (e.g. the cancellation of dynamical pole contributions in the Liischer cor- 
rection method, the general transcendentality structure of the final result) which serve as a 
very strong and nontrivial constraint on the ingredients entering the computation like the 
fine structure of the S-matrix and its most complicated component - the so-called 'dressing 
phase' etc. 

The anomalous dimension of the Konishi operator up to three loops is beautifully 
described by the Asymptotic Bethe Ansatz (ABA) of Beisert and Staudacher within the 
spin-chain approach [1]. 

The four-loop anomalous dimension of the Konishi operator was particularly impor- 
tant as it was only starting from this order of perturbation theory that deviations from 
the ABA appear due to a new topological class of Feynman diagrams - so-called 'wrapping 
interactions'. These effects cannot be understood within the spin-chain point of view and 
their description neccessitates the use of the string a- model in AdS§ x S 5 . The 'wrap- 
ping corrections' are then naturally identified [2] with finite-size Liischer virtual corrections 
charactersitic of a two-dimensional worldsheet QFT. The agreement of direct gauge theory 
perturbative computations [3, 4] and the string cr-model worldsheet QFT computation [5] 
showed effectively that the gauge theory spin chain breaks down exactly in the way expected 
from a two-dimensional QFT on a cylinder i.e. string theory. 

The five- loop order was computed initially only from the string Liischer corrections [6] . 
The main motivation for this computation was that at five loops many more nontrival in- 
gredients of the string cr-model computation start to appear, in particular an infinite set 
of coefficients of the dressing phase and virtual modifications of the quantization condi- 
tions for the constituents of the Konishi state. Here the test of these ingredients were 
initially just the internal cross-checks of the Liischer computation like the cancellation of 
so-called dynamical pole contributions and a simple transcendental structure of the final 
result. A second motivation for this computation was to provide a cross-check for the ex- 
act sets of Thermodynamic Bethe Ansatz equations proposed for the spectrum of J\f = 4 
SYM theory [8-10]. Subsequently full agreement was found, first numerically [11] and then 
analytically [12]. 

It is worth mentioning that the Konishi operator is the first member of an important 
family of so-called twist-two operators, whose anomalous dimensions obey the still myste- 
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rious maximal transcendentality principle [13] and intricate relations with BFKL. In [14] it 
was found that the Asymptotic Bethe Ansatz answer is inconsistent with BFKL starting 
from four loop level, with the discrepancy being attributed to the then unknown wrapping 
corrections. The wrapping corrections for these operators have been computed on the string 
side using Liischer corrections at four- and five-loop level [15, 16], restoring consistency with 
BFKL expectations. 

Liischer corrections for multiparticle states include only the single wrapping contri- 
bution. However the sets of TBA equations [8-10] proposed later, provide a method of 
obtaining the exact anomalous dimension for any coupling thus resumming all multiple 
wrapping effects at once. Numerical computations of the Konishi anomalous dimension for 
a wide range of couplings have been performed in [17] and confirmed in [18]. Subsequently 
large coupling asymptotic fits to the numerical data have been shown to agree with string 
theory computations [19-21]. 

The reason for returning to the six- and seven-loop Konishi anomalous dimension using 
the Liischer correction approach in the present paper has been twofold. Firstly, on the 
gauge theory side new techniques have been developed which allowed for a direct gauge 
theory computation at the five-loop level [22] and the verification of the previous string 
computation [6]. The authors informed us that a six-loop gauge theory computation using 
their methods was within reach and expressed interest in having a string theory cross-check. 
However, numerical methods of solving TBA equations are not precise enough to extract 
such high loop level coefficients, thus making the application of the conventional Liischer 
computation an attractive alternative when applicable. 

Secondly, the infinite set of TBA equations has been recast in a form of a finite set 
of nonlinear integral equations (FiNLIE) [23, 24], which makes a weak coupling expansion 
practically possible - however it is still very complex. This endeavor is now almost finalized 
[25], and the present Liischer computation may provide a nontrivial cross-check. 

For the above reasons we have embarked on computing the six and seven loop anomalous 
dimension of the Konishi operator using the multiparticle Liischer correction formulas of 
[5, 6]. The seven loop level is the theoretical limit of this approach as at eight loops double- 
wrapping will appear which is not taken into account by these formulas. 

The general structure of the result of the multiparticle Liischer correction formula is 



where the integrand can be obtained in a relatively straightforward way. At n-loop the inte- 
grand contains the (n — A) th power of the polygamma function originating from the dressing 
phase. Evaluating the integral by residues [6, 7] derivatives of the polygamma function at 
integers values appear which, after summation for Q, can be expressed in terms of multiple 
zeta values (MZVs). At n-loop MZVs of depth n — 3 appear [7]. Interestingly, irreducible 
MZVs showing up at various stages of the calculations cancel completely leaving only prod- 
ucts of simple zeta functions. We do not have an understanding of this phenomenon yet, 
but will assume that similar structure will appear for the Konishi operator, too. Of course, 
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the fact that our final results can be fitted to a very high accuracy without using irreducible 
MZV's strongly supports this assumption. 

In the present case we have decided to bypass the technical steps above, as at six 
and seven loops they would be prohibitively complicated. We derive just the analytical 
expression for the integrand f(q,Q) and perform the resulting integral and sum semi- 
numerically to a very high precision 1 which allows us to determine the coefficients of the 
zeta-functions in the expected form of the result 

AE^ loop = Cl + c 2 C(3) + c 3 C(5) + c 4 C(7) + c 5 C(9) + c 6 C 2 (3) + c 7 C(3)C(5) (1.2) 

and 

AE 7 _ loop = Cl + c 2 C(3) + c 3 C(5) + c 4 C(7) + c 5 C(9) + c 6 C(H) + c 7 ( 2 (3) + 

cs C(3)C(5) + c 9 C(3)C(7) + cio C 3 (3) + c n C 2 (5) (1.3) 

assuming, as was the case for lower orders, that the coefficients are integers. 

The plan of the paper is as follows. The anomalous dimension of the Konishi operator 
corresponds to the finite volume energy of a two particle scattering state. In section 2 
we derive the ABA momentum and energy of this state up to seven loops. In Section 3 
we express the leading vacuum polarization effects in terms of the scattering matrix and 
evaluate their weak coupling expansion, thus obtaining the integrand of the multiparticle 
Liischer correction formula. These expressions are subject to a consistency check in sec- 
tion 4, where we also compute analytically the highest transcendental part of the answer 
which will be an important cross-check for our numerical approach. In section 5 we discuss 
the semi-numerical computation of our result. Finally we provide the full dimension of the 
Konishi operator up to seven loops in section 6 and make some concluding remarks. 

2 Asymptotic Bethe Ansatz up to 7-loops 

Finite volume energy corrections of multiparticle states have two sources: momentum quan- 
tization and vacuum polarization. Asymptotic BA determines momentum quantization and 
captures all finite size corrections, which are polynomial in the inverse of the volume. In 
contrast, vacuum polarization effects are exponentially supressed and their leading term is 
described by Liischer-type formulas. 

The Konishi operator corresponds to a two-particle state in volume L = 2. These 
particles scatter diagonally on each other and the asymptotical BA equation quantizes 
their momenta as: 

e* L = S(p,-p), (2.1) 

Once the momenta are determined the asymptotic energy of the state can be calculated 
from the dispersion relation 2 

e(p) = ^jl + l&g 2 sin 2 | (2.2) 

1 This is akin to the method used in [16]. In our case, however, a purely numerical summation does not 
suffice to attain the required accuracy and we use several orders of additional asymptotic analysis. 
2 The 't Hooft coupling is reletad to g as g — 
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as 

Eaba(p) = c(p) + e(-p) = 2e(p). 
The scattering matrix takes the form: 

S(p, -P) = <P\ COt l +i e -^ e ^-P) , (2.3) 
e(p) cot | — ? 

where the dressing phase, 6{p, —p), is parameterized in the following manner: 

9{xi,x 2 ) = x{xf, x%) + \{xi,x^) - x(xf,x^) - x(xi,x%) (2.4) 

with 

x ± (p) = -(cot|±i)(l + e(p)) (2.5) 
The functions x have a double expansion 

x(Xl ' X2) = -£^(r-i)t-l) 

r=2 s>r x Jy ' 



_r-l a-1 s-1 r-1 



(2.6) 



and the cv s coefficients, 

CrM = (r-l)(s- l)2cos(^( S - r - 1)) ^ ^zlMj-^M , (2.7) 

are written such that in the weak coupling limit (<? — > 0) the integrals can be easily evaluated 
to yield: 

C2,M = 4Cs5 3 - 40C55 5 + 420Cr5 7 - 4704C 9 <? 9 + • • • (2.8) 

C2, 5 (5) = -8C 5 <7 5 + I68C7/ + • • • 

c 3 , 4 (5) = 24Cs5 5 - 420C75 7 + • • • 
We solve the asymptotic BA equation (2.1) iteratively: 

p = + y - vV + ^V% 4 - 24^3(1 + Cs)5 6 + ^(671 + 960(C 3 + Cs))/ (2.9) 

- — (1269 + I88OC3 + 2O6OC5 + 2100C 7 )5 10 
5 

+ 6^(2643 + 12C 3 (311 + 8C5) + 4072C 5 + 4452C 7 + 4704C 9 )<? 12 

By plugging back this asymptotic momentum into the dispersion relation we obtain the 
asymptotic energy of the Konishi state: 

E BY = 2e(p) = 2 + 12g 2 - 48g A + 336c/ 6 - 12(235 + 24C 3 )g 8 (2.10) 
+12(2209 + 360C 3 + 240( 5 )g 10 
-12(22429 + 4608C 3 + 3672C 5 + 2520( 7 )g 12 

+24(119885 + 24C 3 (1211 + 6C 3 ) + 24156C 5 + 19656C 7 + 14112C 9 )c/ 14 

This provides the exact answer up to g G (three loops) as from four loops onwards wrapping 
corrections start to contribute. 



- 5 - 



3 Liischer type corrections 



The vacuum polarization effects show up in two different ways. First, particles polarize the 
vacuum in which they are moving, this is what we call the direct energy correction. Second, 
the polarized vacuum reacts back and modifies - via the asymptotic Bethe Ansatz - the 
particles momenta, this we call the BA modification. Let us analyze them one after one 
the other. 

3.1 Direct energy correction 

The polarized vacuum will contribute to the energy as 

AE = ~ E / |W^i(9.JP)Sgi(9, -p))e~^ iq)L + 0(g ie ) (3.1) 
Q=l J 

Here we sum up for all bound-states of charge Q of the mirror model, whose momenta are q. 
Sq^ describes how these mirror bound-states scatter on the fundamental "Konishi" particle, 
while eQ(q) denotes their mirror energy. This expression is exact in the weak coupling limit 
up to the order g 14 since at the order g ie double wrapping effects will contribute, too. In 
order to obtain the six- and seven-loop anomalous dimension of the Konishi operator we 
need to expand this expression up to g 14 . Let us analyze the g dependence of the various 
terms. 

Exponential factor 

The mirror energy has the following parametrization 



^ = ^M! ; z± (q,Q) 




± 1 (3.2) 



z+(q,Q) ' Vf "*' 4g \\ qZ + Q* 

Matrix part of the scattering matrix 

There are many representative of the Konishi operator and in the present paper we choose 
the su(2) one which we label by (ll) . This choice has the advantage that only diagonal 
matrix elements will contribute to eq. (3.1). The scattering matrix can be factorized as 



S 



tf$>( Z± > X± ^ = Sscai{z ± ,x ± )S2\{z ± ,x ± )sf i {z ± ,x ± ) (3.3) 

where (a, a) labels the polarization of the mirror bound-states. For one single su(2|2) alge- 
bra the basis for boundstates can be represented in the superspace formalism as (w J 3 w^~ \ 
uP 3 w®~ 2 ~^ 6162, w^w® -1 ^ '6%, w J 3 w^~ 1 ~ :> '62) , where j takes values from to (Q,Q — 2,Q — 
1,Q — 1), respectively. See [5, 26] for the details. As the two su(2|2) factors are equivalent 
the contribution of the scattering matrix can be factored as 

sTr(s2l(q,p)S2l(q, -p)) = S scal {q,p)S sca i(q, -p)sTr(S™$ {q,p)S s rr ^ ) (q, -p)) 2 (3.4) 

We use formulas from [5] and the following identities: 

SZl^x*) = SlZ(x*,z*) = SH^,^) (3.5) 
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where a is the image of a under the transformation 1 -H- 3. 2 f-> 4, but it will not be 
important here as in the supertrace we sum over all states. From [5] we can extract the 
matrix elements, which turn out to be independent of j: 



L„Q-3 



H! W. 



l w 2 



a^(x, z) = Sl(z, x^) 



x^ rji 



x -1- 771 



IK H) 



l w 2 



"2-Ji 



' 3 4 -> 2a|(x,2) = 52(s,x ± ) 



x (1 — x + z ) x 771 / 772 

X + (1 — X~2! _ ) X+ 771 \ 772 



r2 



'3 ->■ al(x, z) = S3(z,x ) 



x 771 772 



771 772 



(3.6) 

(3.7) 
(3.8) 



w\ w, 



h -> i(o§(x,*) + ag(s 1 *)) = 54(^,x±) = s+O glftfc (3 . 9) 

2 z — x^ (1 — x z ) x^ 771 772 



l w 2 



The string frame factors appearing here can be written as ^7 = y f+ an d = jjp- As 

in the Konishi state we have two particles with opposite momenta the factors ^| will cancel: 

sTr(Ci t 2) (^^)Ci?(^-^)) = (0 + 1)%^-^) (3.10) 
+{Q - l)S2(z, x ± )S2(z, -x T ) - QS3(z, x ± ),S3(z, -x T ) - QSA(z, x ± )54(z, -x T ) 
_ 2(x~ - a+)(z+(-l - 2Qx+z~ + z~z+) + x~z-(2Qz+ + x+(-l + z~z + )))z- 
(x" + z~)(z- - x+)(-l + x~z-)(l + x+z-)z+ 

Scalar factor 

The scalar factor for the fundametal particles is 

+ 1 1 

z — x T 1 — ,+„,_ 



X 1 



1 



(3.11) 



and this is what we should fuse for the bound-state: 

S SC al{z,x) = W_S® cal (Zi,: 



(3.12) 



There is a nicer form, however, which does not depend on the individual constituents: 



S- l al (z,x) = Zl 1 (z,x)S su W(z,x) 
The su(2) scalar factor is defined as: 



(z- - x+)(z~ - x-) (1 - ^)(1 - 
while for the case |x | > 1 following [27] we can write 3 

- z log £(3,1(2:, x) = <fj(z + ,x + ) - $(z + ,x~) - <f>(z~,x + ) + $(z~,x~) 

-- [-^(z + ,x + ) + ^(z + ,x~) - ^(z~,x + ) + ^(z~,x~)] 



(3.13) 



(3.14) 



(3.15) 



x-)(x+ - -4^ 2 



(z--x+)(x+-i)(x--i) 



3 We have to be careful as the conventions of [27] are different from ours. To turn into our conventions 



we made the replacements: z 



and x ± -> i^. 
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The functions $ and^ are defined in terms of the integrals 

, . / dw\ f dw 2 1 1 . T(\ + ig(w\ + w7 l - w 2 - w^ 1 )) 
$(x 1 ,x 2 ) = zd ) — — -log 



Cl 2-Ki J Cl 2-Ki wi - x\ w 2 - x 2 r(l - ig{w\ + w 1 1 - wo - w 2 x )) 

(3.16) 

f dw 2 1 , r(l + iq(xi + x-, 1 — wo — Wn x )) 

9(x 1 ,x 2 )=i(b 77-^ log^ ^ S 2" 3.17 

T Cl 2m w 2 -x 2 r(l - ig{xi + x^ 1 - w 2 - w^ 1 )) 

We separate the rational part from the ^ and functions as 

S.cdCs.sc) = S (z,x)e i ^ (3.18) 

where 

„ , , (z~ — x+) 2 (l — z~x~) 

S °^ I » = (^-,-)(, + -x + )(l- a -x + ) <3J9) 

and 

cr(z,x) = 2($(z+,x + ) - 3>(z+,x~) - $(z~,x + ) + $(z - ,aT)) (3.20) 
-(-^(z + ,x + ) + ^(z + ,x~) - ^(z",x + ) + ^(z~,x - )) 

Weak coupling expansion 

Our aim is to calculate the weak coupling expansion of AE for L = 2. In so doing we 
decompose the integrand of the Liischer correction (3.1) 

AE = ~fl [ ^P(q,Qmq,Q) (3-21) 

Q=l J 7F 

into a simpler rational part 

4(x" -x+) 2 (z+(-l - 2Qx+z- + z-z+) + x~z~(2Qz + + x+(-l + z~z+))) 2 fz~ s 



((z+) 2 - (x + ) 2 )((z+) 2 - (x-) 2 )(-l + (X-) 2 (Z-)2)(1 + (X+) 2 (Z-) 2 ) \Z~ 

(3.22) 

which contains both the matrix part and the rational part of the scalar factor, and into the 
more complicated £ part: 

ilog£(<z,Q) = -2($(z+,x + ) - 3>(z + ,x~) - $(z~, x + ) + $(z~, x~)) (3.23) 
-2($(z + , -x") - <$>(z + , -x + ) - $(z", -x~) + $(z~, -x + )) 
-^(z + ,x + ) + ^(z + ,x") - "*(z",x + ) + ^f(z~,x~) 
-^>(z + , -x") + &(z + , -x + ) - &(z~, -x~) + &(z~, -x + ) 

We expand these functions in g 2 as 

P(q, Q) = Ps(q, Q)g 8 + Pio(q, Q)g 10 + Piste, Q)g 12 + Pu(q, Q)g u + ■■■ (3.24) 
£(<?, Q) = 1 + S 2 (g, Q) 5 2 + E 4 (g, Q)s 4 + Se(?, Q)# 6 + • • • (3-25) 
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The expansion of the rational part is quite straightforward. In expanding the ^ and $ 
functions we use the same method we used in [6]. The expansion of the VP(xi, x 2 ) functions 
for |x 2 | > 1 (string region) reads as follows 

*(xi,x 2 ) = - — (^(1 - ig{x! + x^ 1 )) + V(l + ig( Xl + x^ 1 ))) (3.26) 

-^2 (*i(l - ig(xi + xr 1 )) - *i(l + ig(x t + X7 1 ))) 
zx 2 

+9 3 + ^3) (*2(1 - ig(xi + x^ 1 )) + tf 2 (l + ig(xi + x^ 1 ))) 

■ 4 

+^2 ($3(1 - ig( Xl + X7 1 )) - * 3 (1 + i 5 (xi + X7 1 ))) 

DX 2 

—7^ — (^4(1 - ig(xi + xr 1 )) + ^4(1 + igixx + xr 1 ))) + . . . 

where ^ n (x) = (^|) n+1 (log(r(x)) are the standard polygamma functions. If |xi| > 1 then 
$(xi,X2) starts at g e and coincides with ^(xi, x 2 ). In the opposite case using the identity 
$(xi,X2) = 3>(0, X2) — ^(xj" 1 , X2), being valid if |xi| 7^ 1, we can calculate the leading 
expansion of <& as 

9(0, x) = *™ - 2 C39 3 (- + ±) + 24C 55 5 (±+ *)+... (3.27) 

x \x 3x c> / \3x ox y 

$(x7\ x 2 ) = -2C35 3 f % ~ ^) + ■ ■ ■ (3.28) 

V^2 X 2j 

The ^ terms simplify, too: 
ilogE(g,Q)= (3.29) 
2flf(-r - — ) (47^ + #(1 - i#u + ) + ^(1 + 2£ru + ) + ^(1 - igu~) + ^(1 + igu 



X" 1 " X" 

1 



+ * 2 (1 - + *2(1 + ^ + ) + *2(1 - igu~) + * 2 (1 + igu 



o 5 1 1 

^(— - — ) ( 960 (5 + *4(1 - igu + ) + * 4 (1 + igu + ) + f 4 (l - igu~) + * 4 (1 + i^u") 



6 x 
where 

u ± = ^ + (3.30) 

We have to plug back the parametrization of x^ and z and expand in g to obtain £ 2 , 
£ 4 and Eg. This leads to an explicit, however very complicated expression for the direct 
energy integrand in the form 

Pi2(q, Q) + Pio(q, Q)£ 2 (<7, Q) + Ps(q, Q)S 4 (g, Q) (3.31) 
for the 6-loop case and 

Pu(q, Q) + Pi 2 (?, Q)S 2 (?, Q) + Pw(q, Q)S 4 (g, Q) + P 8 (<Z, Q)E 6 (<7, Q) (3.32) 
for 7-loop. 
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3.2 Correction to the asymptotical BA 

The back reaction of the polarized vacuum is a change in the asymptotic BA equation: 

pL + Hog S(p, -p)-2vr = $ (3.33) 

where 

$ = " E / ^M(d q S%l( q ,p))s2l( q , -p))e~^ L (3.34) 
Q=l J 

We can use eq. (3.33) to express 5p as 

* = L + i d p lo g S(p,-p) (3 ' 35) 
and with this correction the change in the energy is 

A P E = 2d p E(p ABA )5p (3.36) 

As the derivative of the dispersion relaton is of order g 2 the momentum shift shows up one 
order higher in the energy formula. Thus it is sufficient to expand $ up to the order g 12 . 
We evaluated the differentiation in the integrand as follows: 

sTr((a g 5g 1 1 (g,p))5g 1 1 (g, -p)) = (3.37) 
(dqS sca i{q,p))S sca i{q, -p)sTt(S s ^ ] {q,p)S^ ] (q, -p)) 2 

+2S sca i{q,p)S sca i(q, -p)sTk({d q S^f (q,p))S s 2af {q, ~p)) sTr (S™£ > (q,p)S^ (q, -p)) 
where 

sTr^Xl'W^Ci?^, = (3.38) 

(Q + l){d g Sl(z, x ± ))51(z, -x*) + {Q- l)(d q S2(z, x ± ))S2(z, -x T ) 
-Q(d q S3(z, x ± ))S3(z, -x T ) - Q(d q SA(z, x ± ))S4(z, -x T ) 

In differentiating the scalar factor we need to differentiate the So part and the a part 
separately: 

d q S scal (q,p) = OVo^x^K^) + S (z,x ± )e i ^ x± Hid q a(z,x ± )) (3.39) 

Weak coupling expansion 

We expand the modification of the BA as 

$ = g 8 $8 + g 10 $ w + g 12 $i2 + • • • (3.40) 
which we write into the form 

$ = " E / Y {dP{q ' Q) + P{q ' Q)dnqi Q)) (3 ' 41) 

Q=l J 7F 
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Besides notations used before we introduced 



dP(q,Q) = (d q So(q,p))S (q,-p)8Ti(S^\z,x ± )S^\z-x^)) 2 + (3-42) 
2S ( q ,p))S ( q ,-p)sTr((d q S^\z,x ± ))S s ^\z,-x^))x 

and 

dS(q, Q) = i(d q a(z, x ± ))S(g, Q) (3.43) 
We need to expand each function to third nontrivial order 

dP(q, Q) = dP(q, Q) 8 g 8 + dP{q, Q W° + dP(q, Q) 12 g 12 (3.44) 

In the a case we need to keep only the first non-trivial term which survives after sym- 
metrization: 

d q <r(z,a*) = ^(-4^-^^)5 g ( 1 I'i(l-^ + )-^i(l+^ + )+ , I'i(l-^")-^i(l+^")) 
2 (x + ) z [x y 

(3.45) 



4 Cross-checks for the integrand 

In the previous sections we have derived the integrand of the six and seven loop contribution 
to the wrapping correction. This integrand consists of two quite distinct parts - the direct 
contribution to the energy which is essentially just the transfer matrix Edirectisi^ Q) an d the 
wrapping contribution to the BY quantization condition E mo d(q,Q). This term cannot be 
expressed in a simple form directly in terms of the transfer matrix and we had to evaluate 
the derivatives of S-matrix elements in order to calculate it (note that in contrast to the 
relativistic case we cannot exchange the derivative w.r.t. the first argument with minus the 
derivative w.r.t. the second argument). We also symmetrized both expressions w.r.t. q. 
The obtained expression has thus the form 

A£>_ ierm = -— ^ / d Q ( E direct(q, Q) ~ E mod (q, Q)) (4.1) 

The above expression is, strictly speaking, just the i^-term Liischer correction, fi-term 
corrections could in principle arise from the residues of the integrand at the dynamical 
poles (s- or t-channel physical poles). However at weak coupling one can argue (see [5]) 
that /i-terms should not appear and hence the residues at the dynamical poles should cancel 
out after summation over Q. This indeed happened in our earlier computations both at 
four and at five loops. This constraint is a nontrivial cross check of the computation of 
both parts of the integrand, since the sum of the residues turns out to vanish only after 
we have included both the direct energy contribution Edi rect (q, Q) and the BY modification 
contribution E mo( [(q,Q). 
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Explicitly we find numerically that 



£ res <?=? 



7» ^direct 

(q,Q) = Yl 2 res g=(? ,£ mod (g,<2) (4.2) 
Q=l g*e{<?5,g6,<?7,<?8} Q =1 9*e{(j5,96,97,g8} 

where the dynamical poles in the upper half plane are 

?5,6 = 5= + i(Q-l) g 7 , 8 = ^i + i(Q + l) (4.3) 

Before we evaluate the full sum and integral numerically in the following section, it is 
also important to derive an independent analytical cross-check for our procedure of numer- 
ically fitting the coefficients of ^-functions for the final expression. 

Experience with earlier four- and five-loop computations show that one can very easily 
derive analytically the coefficient of the single £ function of maximal transcendentality. It 
comes just from evaluating the integral by residue on the 'kinematical pole' q = iQ ignoring 
any other poles coming from the dressing phase. Moreover, this contribution arises only 
from the rational part of the direct energy contribution Edi rec t(q, Q). 

Explicitly, we compute — i res g= jQPi2 (<7, Q), perform a partial fraction expansion and 
look at the term proportional to 1/Q 9 . In the present case we obtain 

_ , 489888 205632 , A A . 

- i res q=iQ P 12 (q, Q) = ^g— + — ^— + • • • (4.4) 

We expect that the coefficient of the maximally transcendental ("(9) is —489888. We will 
not assume it in the subsequent computation but will check at the end whether the numer- 
ically obtained coefficient of £(9) will agree with this predicition. Note, however, that the 
coefficients of £'s with lower transcendentality (like £(7) here) cannot be obtained in this 
way. 

Analogous computations for the seven loop integrand give 

_ , 7318080 
- i res q= i Q Pi A {q, Q) = —qu H • • ■ ( 4 - 5 ) 

thus predicting that the coefficient of C(H) in the final answer will be 7318080. Again, we 
will not use this information in the numerical fits, but will check whether this coefficient 
arises independently from the numerics. For the seven loop integrand we have also checked 
numerically the cancellation of dynamical poles. 

5 The evaluation of the multiparticle Liischer correction 

We will now proceed to evaluate the expression (4.1). We will first give details for the six 
loop computation and then briefly comment on the modifications at seven loops. 

It is easy to see that a natural variable for the integrand is s = q/Q and one can see 
that the integrand approaches at large Q a single scaling function of s (times an appropriate 
power of Q). Hence we introduce 

fs-ioop(s, Q) = Q(E direct (sQ, Q) - E mod (sQ, Q)) (5.1) 
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and evaluate 

ro ° ds 



- Y, / -h-ioop(s,Q) (5.2) 

We first numerically evaluate using Mathematica the integrals from Q = 1 up to Q ma x — 
2000 using 80 digit precision (we set PrecisionGoal to 80 and WorkingPrecision to 100). 
By itself, this is a far too small Q ma x to estimate the infinite sum with sufficient precision in 
order to reliably fit the coefficients of the £'s. We therefore supplement this numerical part 
of the computation by computing analytically the asymptotics of the integral around Q = oo 
to several orders and performing the sums on these asymptotics from Q — Qmax 

+ 1 = 2001 

to Q = oo. 

Explictly we find that the six-loop integral has the following asymptotic expansion at 
large Q: 

f6-loo P (s, Q) ~ Q9T^r (o-n(s) + b n (s) log ^ + c n (s) log 2 (5.3) 

We evaluated the first six orders (up to 1/Q 19 ) and integrated the coefficient functions 
On(s), b n (s), c n (s) (two of these integrals could be done analytically by Mathematica but 
80-digit numerical evaluation was faster and just as reliable). In this way we obtained the 
large Q expansion of the final summands 

U n + B n log^- + C n log 2 ^-) (5.4) 



As a cross check, we verified that evaluating the above expression at our cut-off point 
Q = 2000 agrees with the numerical integral for the same Q up to 10 -61 . Now we sum the 
asymptotic expression (5.4) from Q = 2001 to oo and add it to the numerical sum. 

We use EZ-Face [28] to express the resulting number as a linear combination with 
integer coefficients. Using 44 digits precision we obtained 

A£ 6 _ /oop = 261468 - 207360 C(3) + 156384 C(5) + 105840 C(7) - 489888 C(9) 

+155520 C(3)C(5) - 20736 C 2 (3) (5.5) 

We see that the the coefficient of £(9) obtained from the above semi-numerical procedure 
exactly agrees with the analytical prediction of previous section. 

In fact our numerical expression for AE obtained using Qmax — 2000 and subsequent 
six orders of asymptotic expansion agrees with the above analytical expression up to 10 -59 
(this is more than 60 decimal digits agreement as AE = —42127.1014157...). In order to 
make a further test, we evaluated an additional seventh order of the asymptotic expansion. 
Incorporating its contribution to the numerical expression for AE increased its agreement 
with the analytical expression (5.5) up to 10 -66 . Finally, let us note that if we were to 
keep only the numerical sum up to Q max = 2000 without the asymptotic expansions, the 
precision would be much too small to fit the £ functions. Indeed the deviation of this sum 
from the exact answer starts already at order 10 -20 . 



With hindsight, setting AccuracyGoal would be a more efficient option. 
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The evaluation of the seven loop wrapping correction proceeded along the same lines. 
We again evaluated the integrals numerically from Q = 1 to Q = 2000 using this time 
slightly different settings (we set now AccuracyGoal to 70 and WorkingPrecision to 90). 
We could not increase these settings due to some quirks in Mathematica 8. Indeed, at the 
end of the computation we found that Mathematica 8 did not evaluate the integrals for 24 
values of Q, so we had to evaluate these remaining integrals using Mathematica 7. 

Subsequently, we evaluated 9 orders of asymptotic expansion, the reason being that for 
the seven loop case we needed to fit 11 coefficients in contrast to 7 at the six loop level, 
and hence we needed maximal numerical accuracy. We thus evaluated 



in a similar manner as outlined above. The above expression evaluated for Q = 1999 agreed 
with the numerical integral for that value of Q up to 10 -75 . We then summed the above 
asymptotic expression from Q = 2001 to Q = oo and added the result to the sum of the 
numerical integrals computed earlier. 

In order to fit the coefficients of the £'s, h is advantageous to reduce the magnitude 
of the coefficients by factoring out an integer divisor. Experience with lower loop levels 
suggests to factor out 36. Now 65-digit precision in EZ-Face [28], allowed us to obtain the 
analytical answer: 



A£ 7 _zoo P = 36 • - 139230 + 125904 C(3) - 22472 C(5) - 47964 C(7) - 32616 C(9) + 



The deviation of the above expression from our semi-numerical answer appeared only at 
order 10~ 73 , therefore right at the edge of our numerical precision. However, a key test of 
the above result is the coefficient of £(11) which was fitted to be 36 • 203280 = 7318080 and 
exactly coincides with our analytical prediction from the previous section. 

6 The final result and conclusions 

In order to compute the final anomalous dimension of the Konishi operator, it remains to 
add the six and seven loop wrapping corrections (5.5) and (5.7) to the relevant terms in 
the Asymptotic Bethe Ansatz answer (2.10). For completeness, we will present here the 
full result up to seven loops, incorporating also the four loop wrapping correction obtained 




(5.6) 



+ 203280C(H) - H808C 2 (3) + 11424C(3)C(5) - 53760 C(3)C(7) + 




(5.7) 



in [5] 



AE 4 _ loop = 324 + 864 C(3) - 1440 C(5) 



(6.1) 



and the five loop one from [6] 




11340 + 2592 C(3) - 5184 C 2 (3) - 11520 C(5) + 30240 C(7) 



(6.2) 
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The final answer for the Konishi dimension is 

A = 4 + I2g 2 - 48s 4 + 336g 6 + 96(-26 + 6 C(3) - 15 C(5))s 8 (6.3) 
-96(-158 - 72C(3) + 54C(3) 2 + 90C(5) -315C(7))s 10 

-48(160 + 432 C 2 (3) - 2340 C(5) - 72 C(3)(-76 + 45 C(5)) - 1575 C(7) + 10206 Q{$))g 12 

+48(-44480 -8784C 2 (3) + 2592C 3 (3) - 4776C(5) -20700C 2 (5) 

+24C(3)(4540 + 357C(5) - 1680 C(7)) - 26145 C(7) - 17406 C(9) + 152460 C(ll))# 14 

Our six loop result agrees with the computation of [25] obtained from FiNLIE thus being a 
very nontrivial cross check of the proposed finite set of nonlinear equations describing the 
spectrum of M = 4 SYM. It would be very interesting to compare this result with a direct 
gauge theoretical computations along the lines of [22]. It would be also very interesting to 
analytically verify our seven loop result, as in this case the identification was performed 
at the edge of computational precision (however with a nontrivial consistency check of 
the analytically known coefficient of ("(11)). Finally let us note that extending the result 
for Konishi to eight loop order would be conceptually quite important, as at that order 
double wrapping effects will appear. These phenomena were analyzed for the nonprotected 
vacuum state in the 7-deformed theories in [29] but they are beyond the reach of the Liischer 
correction formulas for excited states known up to now. 
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